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^ . Abstract 

^' 

We calculate the low-energy limit of heavy quark current correlators within an expansion 
in the inverse heavy quark mass. The induced low-energy currents built from the gluon 
fields corresponding to the initial heavy quark currents are obtained from an effective 
action for gauge fields in the one-loop approximation at the leading order of the 1/m 
expansion. Explicit formulae for the low-energy spectra of electromagnetic and tensor 
heavy quark current correlators are given. Consequences of the appearance of a nonvan- 
ishing spectral density below the two-particle threshold for high precision phenomenology 
of heavy quarks are discussed quantitatively. 



1 Introduction 



Quantum corrections can qualitatively change the analytic structure of the Green function 
or the symmetry properties of a classical field-theoretical system. Qualitatively new features 
compared to the tree level picture may emerge already at finite orders of perturbation theory 
for Green functions while some effects can only appear as a result of summing up an infinite 
number of perturbative terms. An effect of non- conservation of the Abelian axial current 
reveals itself at leading order as a result of the calculation of an one-loop triangle diagram |jl], ^ 
while spontaneous symmetry breaking or bound state formation cannot be observed at any 
finite order of perturbation theory and requires an infinite summation of relevant subsets of 
diagrams. For problems related to an investigation of the symmetry properties of quantum 
systems, such an infinite summation can be readily done by introducing an effective action 
for the system and calculating it as a loop expansion series. Such an approach reorders the 
perturbation series with respect to the lines of diagrams related to external fields and allows 
one to take into account the entire dependence on external fields exactly within a given order 
of the loop expansion. An effective action can be considered as a generating functional for the 
vertex (one-particle irreducible, or proper) Green functions (see e.g. Ref. 0). The treatment 
of external fields beyond finite order of expansion was done in Ref. [Q. The efficient method of 
calculating the effective action is based on the Legendre transform of the generating functional 
for connected Green functions and heavily uses the functional techniques . One can also use 
a practical calculation technique by substituting the shifted fields in the original Lagrangian of 
the system |^. The part of the effective action which is constructed from the constant field is 
usually called an effective potential and is used to analyze the fundamental symmetry properties 
of the theory beyond the plain perturbation theory where the effect of the external fields is 
resummed to all orders. The expansion leading to the effective potential is, in fact, an expansion 
in Planck's constant h, i.e. the correction accounts for the deflection from the classical limit. 
Therefore new effects which are absent in the classical approximation may appear within such 
an approach. An example for this kind of new quantum effects is the light-by-light scattering 
which emerges as a quantum correction to the photon dynamics due to the interaction with 
virtual electrons. In the low-energy limit it can be seen as a nonlinearity of the equations for the 
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strong electromagnetic fields in the vacuum. The behaviour of the electromagnetic fields with 
such a correction can be described by the Euler-Heisenberg Lagrangian The generalization 
to non-Abelian fields is discussed in Ref. 0. The effective potential is also a powerful tool for 
investigating effects of spontaneous symmetry breaking by quantum corrections [|10| and for 



analyzing the properties of particle systems at finite temperature and density [11 



A special advantage of using the external field technique in gauge theories is an explicit 
gauge invariance of the effective action that allows to drastically simplify the computation 
and to reduce the number of necessary diagrams [0]. It is generally believed that in non- 
Abelian gauge theories the nonperturbative fluctuations (instantons) create a complex 
vacuum structure that eventually explain (or is responsible for) the low-energy observable 
spectrum of particles . The technique of calculating in external fields was heavily used for the 
calculation of quantum corrections to the effective action of gauge fields in the classical instanton 



background [|T^. In practical applications to QCD the complex vacuum structure could explain 
a phase transition from the quark-gluon representation of Green functions at high energies to the 
hadron picture at low energies. While the problem of a full description of this transition remains 
unresolved, Wilson's operator product expansion which is one of the key tools for calculating 
the correlation functions at short distances is used for describing hadronic properties at low 
energies in a semiphenomenological way using sum rule techniques |T^. The external field 
technique provides a convenient way for practical calculations within the sum rule method 



based on a semiphenomenological account for the condensates of the local operators [|T^, [T^ . 

In the present paper we calculate the low-energy limit of heavy quark current correlators 
within the 1/m expansion where m is a heavy quark mass. The induced low-energy currents 
corresponding to the heavy quark initial operators are obtained from an effective action for 
gauge fields of the SU{Nc) ® U{1) model in the one-loop approximation at leading order of 
both the coupling constant and the 1/m expansion. Our results for the effective action for 
vector and tensor currents are presented in Sec. 2. Using the external field technique as a 
convenient framework for practical calculations of the induced currents in this model, in Sec. 3 
we present these induced currents. In Sec. 4 we show phenomenological applications, and in 
Sec. 5 we deal with the consequences encountered for calculating arbitrary moments. 
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2 The effective action 



While the technique is standard, we are aiming at concrete results for further phenomenological 
applications to sum rules for the vacuum polarization functions of heavy quarks. Therefore we 
briefly outline the calculation and the related issues. Further details can be found in Refs. 
The Lagrangian of a heavy fermion field interacting with a gauge field B of the gauge group 
SU{N^)®U{i) reads 

L = i,{iYd^ + ^^B^-m)ij (1) 

where B^ = eA^ + gsB^. Here is a gauge field of the f/(l) subgroup (photon) with the 
coupling constant e and 5^ is a gauge field of the SU{Nc) subgroup (gluon) with the coupling 
constant g^. The matrix notation for the non-Abelian gauge field potentials is used, = taB"^, 
ta are generators of the gauge group SU{Nc). A generating functional W[J] of connected Green 
functions is given by a functional integral with the sources J, 



Z[J] = exp{m[J]) = j[di)di)]ey.^(^ j L\x)d^x^ 

[# d%l)] exp {i J [ipiiYdf, + l^B^ - m)i) + J^B^ d^x^ (2) 



where the product J^B^ implies a trace with respect to the representation of SU{Nc) U{1). 
A proper gauge fixing is implied as well. The effective action T[B] for the gauge field is then 
given by the Legendre transform 

T[B] = W[J] -JB, B= (3) 

oJ 

It was shown that this procedure is equivalent to the more direct calculation in external 
fields (see e.g. Ref. [0). It is also a generalization of results obtained for constant external 
fields 0. Up to leading order in h the effective action constructed with a Legendre transform 
can also be found through 

exp{{r[B]) = j [# d^l)] exp J ^(^7^9^ + 7^i3^ - m)V d^x^ = det {i^d^ + ^^B^ - m) (4) 

where B is now a classical gauge field (integration over d^x is implied for expressions of the 
action). By using the identity det M = exp(tr(lnM)) for an operator M one continues with 

iV [B] = tr [ In {i^d^, + ^^'B^ - m) 1 . (5) 
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Figure 1: Heavy quark loop correction to the electromagnetic current 



2.1 Results for the effective action 

A straightforward calculation of the functional determinant gives a correction to the effective 
low-energy action of the gauge fields within SU {Nc)^U{l). We calculate the leading nontrivial 
contribution in the 1/m expansion. The functional determinant in loop expansion can be 
represented by Feynman diagrams. A diagram which gives a correction to the effective action 
due to a heavy quark loop is shown in Fig. ||. Two-gluon transitions are forbidden according 
to a generalization of Farry's theorem to non-Abelian theories ||2^. We are interested in the 
behaviour of the amplitude associated with the diagram in Fig. |l| at low energies and, therefore, 
take the limit of a very heavy quark. Formally the limit m — oo is taken which in physical terms 
means that m is much larger than all momenta of external legs of the diagram in Fig. namely 
the three gluons and the photon. A straightforward calculation of the diagram presented in 
Fig. gives the one-loop result for the correction to the effective action induced by a heavy 
quark loop at leading order of the 1/m expansion. It reads 



where F^j, is a field strength tensor for the f/(l) subgroup, Gf^u = taG'^^j is the field strength 
tensor for the SU{Nc) subgroup, and dabc are the totally symmetric SU{Nc) structure constants 
defined by the relation date = 2tr({)f:a, tb}tc)- The trace in Eq. (H) is understood as a trace with 
respect to the Lorentz indices of the fields, i.e. one considers the strength tensors of gauge 
fields as matrices for which tr(FC) = F^^G^^. This makes the formulae shorter and more 
transparent. Before we proceed to the calculation of the induced current in the next section, 
we compare the result with the corresponding expression in QED. 



AF, 



^9s dabc 



-[Uii{FG''G^G^) - 5tr(FG"^)tr(G^G^) 



(6) 



QCD — 



180m4(47r) 
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2.2 Comparison with QED 

The effective action within QED corresponding to Eq. (^) is known as Euler-Heisenberg La- 
grangian p!9| , 



Ar 



2a' 



QED 



{E' - ny + 7{E ■ H) 



a 



Arc 



(7) 



which can be obtained from Eq. (H) by replacing dabc by 1 and G^j, by Ffj_^, supplemented by 
the obvious substitution Qs e. Using the definitions for the electric field E and the magnetic 
field i?, 



-E„ F 



iO 



one finds 



iO 



E- F'^ = F ■ 



tr(F2) = 2(^2 - #2), tr(F^) = 2{E'^ - + A{E ■ Hf, 
ii{FF)=A{E-H), tr(FFFF) = -^(tr(F2))2 + tr(F^) =4(E-i?)2 (9) 



where 



F = -e FP" 



(10) 



Note for further use and convenience that in order to rewrite the expressions from one form to 
another there is one more relation between the fourth-order monomials of the electromagnetic 
field strength tensor F^,^, 



(FFY = -2{Fy + AF' 



with 



FF = F>"'F, 



pA ^ F^^'^F^^F^^Fp^. 



(12) 



Therefore the correspondence between Eq. (^) and its QED counterpart in Eq. (0) is established. 



3 The induced current 



An expression for the induced electromagnetic current as being an effective electromagnetic 
current for the low-energy effective theory describing the interaction of photons and gluons is 
given by the derivative of the effective action with respect to the external Abelian gauge field, 



5T[B\ 
' 5A„ 



-tr 



(13) 



The derivative with respect to can be replaced by a derivative with respect to = 

dA^ Ohf^fyi dhy^ 
3.1 Results for the induced vector current 



With the exphcit expression for the effective action given in Eq. we obtain [21 



= 9^(9^^ C"^^ = — L f [14(G"'G''G")^'' - h{G''Y''ii{G^G% (15) 

Note that the current conservation dyJ^^ = is automatically guaranteed because the operator 
O^^ is antisymmetric, O^'^ + O^^ = 0. Higher order corrections in the coupling constant a of 
the U{1) subgroup are omitted. The induced electromagnetic current in Eq. ( ]r5| ) is a correction 
of order 1 / in the inverse heavy quark mass which vanishes in the limit of an infinitely heavy 
quark. Corrections in the inverse heavy quark masses are important for tests of the standard 
model at the present level of precision and have been already discussed in various areas of 



particle phenomenology 22, 23, 24 . 



3.2 Other types of induced currents 

Expressions for the induced currents with quantum numbers other than that of the electro- 
magnetic current J^*" = 1 can be obtained in a similar way. The axial current of fermions 
naturally appears in the standard model as a result of an axial-vector interaction of fermions 
with W and Z bosons. The corresponding induced current can be obtained as a derivative of 
the respective effective action with respect to the Z boson field. The leading order diagram, 
however, contains only two external legs and is UV divergent. In case of massless fermions 
this diagram leads to the anomalous non-conservation of the axial current that also requires a 
strict definition of the corresponding operator ip'^^'^^ip within perturbation theory because its 
renormalization is not dictated by the Ward identities any more. The leading corrections to the 
anomaly of the axial current for massless fermions due to strong interactions were considered in 
where implications of the Adler-Bardeen no-go theorem (non-renormalizability) 
were discussed also for dimensional regularization and different definitions of the axial current 
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at the tree level. Higher order corrections were analyzed in Ref. [^]. Explicit high order cor- 
rections to the expression for the anomaly depend on the renormalization prescription for the 
composite operator ip'j^'y^ip within perturbation theory. 

The scalar current ipip appears in an interaction vertex for the Higgs boson and was inten- 
sively studied. The decay if — 77 is described by the effective interaction 

ALh = gsaHFF (16) 

with the parameter gs (remark the different to Qs, S stands for the scalar current) given by 
the corresponding one-loop diagram with two external photons. Here H is an interpolating 
field for the Higgs boson. The decay of the Higgs boson into two photons is calculated up to 



high orders of perturbation theory and mass expansion ||28|. The application of the effective 



potential technique to the analysis of the correlators of the scalar gluonic currents which emerge 



in the decay of the Higgs boson into hadrons was considered in Ref. [^]. 



3.3 Results for the induced tensor current 

In the context of the effective action we consider a tensor current of the form 

JA''^ = ^a'^^^, a^"' = ^[Y,Y] (17) 

and calculate its low-energy limit induced by a heavy quark loop. The properties of this current 
are rather similar to those of the electromagnetic current. Note that the classical vector mesons 
[p, uj, 0) interact with this current and can be produced by it. We introduce an interaction 



ALt = gr^a'^'tpF, 



in the Lagrangian of heavy quarks and readily find the effective action for gauge fields induced 
by such a vertex. The low-energy limit at the one-loop order reads 

Tt = ^gQ^f^4^'j^2 (32tr(FG"^G'G'^) - 15tr(FG'')tr(G''G^)) (19) 

with the same notations as in Eq. (^. According to the form of the effective interaction in 
Eq. (|1^) the induced current J^^, is given by a derivative 

StJ'" = -jf- (20) 



and explicitly reads 



J^"" = ^^^s^^ly (32{G^G''G')'"' - 15{GY"^^iG''G")) ■ (21) 
Note the lower power of the heavy quark mass in Eq. (|21|) as compared to Eq. ([I5|). 



4 Phenomenological applications 

High precision tests of the standard model remain one of the main topics of particle phenomenol- 
ogy |3^. The recent observation of a possible signal from the Higgs boson may complete the 
experimentally confirmed list of the standard model particles [0. Because experimental data 



are becoming more and more accurate, the determination of numerical values of the parameters 
of the standard model Lagrangian will require more accurate theoretical formulae. Recently 
an essential development in high-order perturbation theory calculations has been observed. A 
remarkable progress has been made in the heavy quark physics where a number of new physical 
effects have been described theoretically with high precision. The cross section of top-antitop 
production near the threshold has been calculated at the next-to-next-to-leading order of an 
expansion in the strong coupling constant and velocity of a heavy quark with an exact account 
for Coulomb interaction (as a review, see Refs. [^, |3^). This theoretical breakthrough allows 
for the best determination of a numerical value of the top quark mass from experimental data. 
The method of Coulomb resummation resides on a nonrelativistic approximation for the Green 
function of the quark-antiquark system near the threshold and has been successfully used for 
the heavy quark mass determination within sum rule techniques [0, Being applied 

to quarkonium systems this method is considered to give the best estimates of heavy quark 
mass parameters ^ ^ Technically an enhancement of near-threshold contributions 



to sum rules is achieved by considering integrals of the spectral density of the heavy quark 
production with weight functions which suppress the high-energy tail of the spectrum. The 
integrals with weight functions 1/s" for different positive integer n, s = E"^, where E is the 
total energy of the quark-antiquark system, are called moments of the spectral density and are 
most often used in the sum rule analysis The interest in the precision determination of 



the c-quark mass is especially high because this parameter introduces the largest uncertainty 
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to the theoretical calculation of the running electromagnetic coupling constant at Mz which is 



one of the key quantities for the constraints to the Higgs boson mass We stress that in the 
corresponding determination of the running electromagnetic coupling constant at Mz based on 
direct integration of experimental data over the threshold region, the sensitivity of the results 



to the c-quark mass is much weaker [43|. 

For phenomenological applications of our results one therefore has to compute the two- 
point correlation functions of induced currents. In the following we show that there is a strong 
constraint on the order n of the moment that can be used in heavy quark sum rules. Because 
of the contribution of low-energy gluons, only the first few moments exist if the theoretical 
expressions for the correlators include the order of perturbation theory. 



4.1 The spectrum for the induced vector current correlator 

First we discuss the case of the vector current where the data are obtained from e^e~ annihila- 
tion experiments and are rather precise. The basic quantity for the analysis of a vector current 
= g7^g of a heavy fermion q within sum rules is the vacuum polarization function 

127r2^ j {TUx)MO))e'''d'x = {q,q, - g,.q^mq^). (22) 

With the spectral density p(s) defined by the relation 

p(s) = -^(n(s + zO) -n(s-zO)), s>0, (23) 

the dispersion representation 

n(,^)=/^^^^^ (24) 

J s — q"^ 

holds. A necessary regularization and subtraction is assumed in Eq. (|^). The normalization 
of the vacuum polarization function 11 (g^) in Eq. ( p^ ) is chosen such that one obtains the high- 
energy limit lims_^oo p{s) = 1 for a lepton. For the quark in the fundamental representation of 
the gauge group SU{Nc) the high-energy limit of the spectral density reads p{oo) = Nc. The 
integral in Eq. (|2^ runs over the whole spectrum of the correlator in Eq. ( p^ or over the whole 
support of the spectral density p{s) in Eq. (pSl). 

A correlator of the induced vector current has the general form 

{TJ^'{x)r{0)) = -d^dfs{TO^"'{x)0'''^) (25) 
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Figure 2: Induced massless correlator diagrams 



where an explicit expression of the current as a derivative of the antisymmetric operator O^"^ 
has been employed. The resulting correlator {TO^°'{x)0^^) in Eq. (25) contains only gluonic 



operators. Such correlators were considered previously in the framework of perturbation the- 
ory p4| , |45[] . In leading order of perturbation theory the correlator in Eq. ( pS] ) has a topological 



structure of a sunset diagram, as it is shown in Fig. 0(a). Technically, a convenient procedure 
of computing the sunset-type diagrams is to work in configuration space P6|. We find 



A Fourier transform of the correlator in Eq. ( 126|) gives the vacuum polarization function in 
momentum space which reads 

12txH I {TJ^ix).MO))e''^d''x = (q^q, - g.^q'^Mq^) (27) 

where at small q^ {q^ <^ rn?) 

For QCD with the colour group SU{?>) one has dabcdabc = 40/3. The spectral density of the 
vacuum polarization function n(g^) in Eq. (pTl) is given at small values for s by 

P(^)U = C7,(£^)'. (29) 



Note that the spectral density in Eq. ( PUj ) can be found without an explicit calculation of its 
Fourier transform. Instead one can use a spectral decomposition (dispersion representation) in 
configuration space which was heavily employed for the analysis of sunset diagrams in Ref. ||46|| . 
In this particular instance the spectral representation of the correlator in configuration space 
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reads 

i vr^ 



28r(6)r(5) Jo 

with Di^x^, s) being the propagator of a scalar particle of mass -y/i 



/ s^Dix^,s)ds (30) 
Jo 



where Ki{z) is the McDonald function (a modified Bessel function of the third kind, see e.g. 
Ref. [0). r(z) is Euler's gamma function. 

An asymptotic behaviour of the spectral density of the corresponding contribution for large 
energies (where the limit of massless quarks can be used) enters the expression for the ratio 
R{s) of e^e~ annihilation into hadrons and has been known since long ago [|^, ^ |50[. This 
term is usually called light-by-light (Ibl) contribution and reads 

^iM(^) = )' ^ - 128C(3)) . (32) 

Here C,{z) is the Riemann C, function with (^(3) = 1.20206. . . The contribution to the spectral 
density given in Eq. (|3^) is negative while our result given in Eq. (^9]) is positive as it should 
be the case for the spectral density of the electromagnetic current as a Hermitean operator. 

4.2 The spectrum for the induced tensor current correlator 



The results for the correlator of the tensor current given in Eq. (|2l|) are slightly more compli- 
cated. The correlator reads 



U-K^i J {Tjf,^{x)jc,f3{0))e'''''d^x = {gf,a9ui3 - 9,ip9ua)'n.g{q'^) 

+i9,iaquqf3 - 9t^(sq^qa - 9uaq,iqt3 + 9^fsqtMqa)'^qiq'^) (33) 



where two scalar amplitudes are possible now. One finds 



q^ , 2n Qldabcdabc / q^ V, / /i^ 



The physical content of the amplitudes Uglq"^) and nq(g^) is related to contributions of the 
states with J^*^ = 1 and J^^ = 1^ , resp. Note that the sum rule analysis for the mesons 
with quantum numbers J^*" = 1^ has been done in Ref. |51| with quark interpolating cur- 



rents. From the present results we also see a possibility to use gluonic currents as interpolating 
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operators for such mesons. The vahdity of such a description depends strongly on the strength 
of the interaction of the meson in question with the corresponding interpolating operator which 
is difficult to estimate independently. 

Note that there are only two independent gluonic operators to construct the induced currents 
under consideration. The electromagnetic current is given by a derivative of a special linear 
combination of these operators while the tensor current is given by a linear combination of 
operators themselves. There is one more current relevant to the situation. It originates from 
the Gordon decomposition of the electromagnetic current (see e.g. Ref. ||l9| ) 



2m7/'7^V' = d^i'ipa'^'ip) + tp iD^'ip. (35) 

This relation retains for the induced currents as well. The left hand side and the right hand 
side of Eq. ( |5BD have different parity as for the number of Dirac 7-matrices between spinor fields 
which is refiected in an additional factor m at the left hand side of Eq. (^). In the massless 
limit these types of currents are alien and can never mix. At the level of induced currents the 
Dirac structure of the initial heavy quark currents is refiected in different degrees of suppression 
by the heavy quark mass m. 

4.3 The spectrum for a mixed current correlator 

Having both currents at hand, one can study a mixed correlator of the form 

Un^t I {TJ^{x)J^f,{0))e''^d^x = i{g^^qp - g^pq^)IiM{q^) (36) 
with a single scalar amplitude nj\,/(g^). Such mixed correlators are useful in sum rule applica- 



tions |£2|. One finds 



2x Q7dabcdabc fas\^-q / 9 \ i„ / 



The physical content of the amplitude is given by the J^*^ = 1 resonance, i.e. by the T-family 
in case of b quarks. 



13 



5 Moments of the spectral density 

As mentioned earlier, the moments of the spectral density p(s) of the form 



M.-r-^ (38) 



are usually studied within the sum rule method for heavy quarks These moments are 

related to the derivatives of the vacuum polarization function 11 (g^) at the origin, 



(39) 



Such moments are chosen in order to suppress the high energy part of the spectral density p{s) 
which is not measured accurately in the experiment. Within the sum rule method one assumes 
that the moments in Eq. (|38D can be calculated for any n or, equivalently, that the derivatives in 
Eq. ([39| ) exist for any n. The existence of moments seems to be obvious because one implicitly 
assumes that the spectral density p{s) of the heavy quark electromagnetic currents vanishes 
below the two-particle threshold s = 4m^ which means that the vacuum polarization function 
n(g^) of heavy quarks is analytic in the whole complex plane of except for the cut along 
the positive real axis starting from 4m^. This assumption about the analytic properties of 
the vacuum polarization function n(g^) is known to be wrong if a resummation of Coulomb 
effects to all orders of perturbation theory is performed: as a result of such a resummation the 
Coulomb bound states appear below the perturbation theory threshold s = Aw?. 

5.1 Infrared singular behaviour of the moments 

The qualitatively new feature of effective currents given in Eqs. (|T5|) and (^l]) is that they are 



expressed through massless fields. Therefore the spectrum of the two-point correlators of these 
currents start at zero energy. This feature drastically changes the analytic structure of the two- 
point correlators of these currents and, in particular, their infrared (IR) or small s behaviour 
because of the branching point (cut) singularity of n(g^) at the origin = 0. This new 
feature of having a nonvanishing spectrum below the formal tree-level two-particle threshold 
which appears at the 0{a^) order of perturbation theory for induced current correlators has 
important phenomenological consequences. Indeed, such a change of the analytic structure 
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of induced current correlators affects strongly the theoretical expressions for some observables 
usually employed in heavy quark physics for the precision determination of the parameters of 
heavy quarks and their interactions. 

Because of the low-energy gluon contributions, the large n moments of the spectral density 



in Eq. (|38|) do not exist and cannot be used for phenomenological analyses. Caused by the factor 
{s/Am?Y Eq. (^) the moments become IR singular for n > 4 in case of the induced vector 
current. This can already be seen by looking at the factor I/tjv" in the induced vector current 
in Eq. (|15|) . For the induced tensor current the corresponding moments start to diverge earlier 
because of a weaker suppression by the heavy quark mass, the corresponding factor in Eq. ( pi| ) 
is instead of Therefore, in this case the moments become IR singular already for 

n > 3. Note that in early considerations of sum rules quite large n were used. For instance, the 
numerical value of the gluon condensate was extracted from sum rules for the moments with 



n ~ 10-7-20 In view of our result on the low-energy behaviour of the spectral density, 

one has either to limit the accuracy of theoretical calculations for the moments to the 0{af) 
order of perturbation theory which seems insufficient for a high precision analysis of quarkonium 
systems (especially if the Coulomb resummation to all orders is performed) or to use only a 
few first moments. For small n, however, the high-energy contribution, which is not known 
experimentally with a reasonable precision, is not sufficiently suppressed and introduces a large 
quantitative uncertainty into the sum rules for the moments. An analysis based on finite energy 
sum rules is free from such a problem and can be used in phenomenological applications ||54 . 



Note in passing that there is no low-energy gluon contribution (and therefore no low-energy 
divergence problem) for correlators of the currents containing only one heavy quark with mass 
m. The spectrum of such correlators starts at m? and there are no massless intermediate 
states contributing to the correlator in perturbation theory (see e.g. Ref. [^). The theoretical 
expressions for such correlators can be used for high precision tests of theoretical predictions 
when the accuracy of experimental data in corresponding channels will improve in the future. 
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5.2 An infrared safe moment definition 

The infinite-integration sum rules with large n can be retained at high orders of perturbation 
theory if an appropriate cutoff at small energies is introduced. This can be readily achieved by 
calculating the moments (or derivatives) at some Euclidean point = —A ||56[. Indeed, for 
the regularized moments 



(40) 



g2=-A J (s + A)"+i 

there is no divergence at small s. However, the regularization parameter A cannot be arbitrary 
small. The reason is that the resulting correlator of gluonic currents in Eq. ( P^j ) is essentially 
normalized at /i^ = A when radiative corrections are taken into account, like the diagram shown 
in Fig. 0(b). The corrections are known to be large and, therefore, A is much larger than the 
expected IR scale in quark channels [^. This observation makes the phenomenological analysis 
based on the regularized sum rules in Eq. (^) unprecise even for reasonably large n because the 
continuum contribution to moments is not suppressed for large values of A. The suppression 
of the high-energy tail can be enhanced by constructing a more special kind of moments which 
exploit the explicit behaviour of the spectral density at low energy. The IR safe moments with 
a more efficient suppression of high-energy tail for the induced vector current have the form 

It seems that even in such an optimized form these moments cannot be used for the cc-system, 
i.e. for analyzing J/ip resonances because the suppression is insufficient for achieving the preci- 
sion goals necessary for the c-quark mass determination. For the 66-system, i.e. for the analysis 
of T-resonances the dependence on the cutoff A is numerically essential for moments at large 
n and for A as large as A = 1 2 GeV^ which is still small according to the estimates of 
radiative corrections in gluonic channels. We would like to remind the reader that Coulomb 
poles which are essential for the analysis of the T-resonance and 66-production near threshold 
give contributions which are formally of the 0{al) order (for the value of the Coulombic wave 
function at the origin see e.g. Ref. |3^), an order which coincides with the order of corrections 
considered here (see Eq. (^8])). The theoretical expressions for the correlators in the scalar 
channel where below threshold corrections start at the 0(q;^) order are more sensitive to these 
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special below threshold contributions. However, data in the scalar channel are considerably 
worse than those in the vector channel and there is no possibility of a high precision analysis 
on the scalar channel at present. 

6 Conclusions 

We have presented corrections to the heavy quark currents induced by a virtual heavy quark 
loop and expressed through the gluon operators. We have considered the electromagnetic 
current and the tensor current closely related to it. Heavy quark loop induced corrections first 
appear at 0{ag) order of perturbation theory and are given by the term of the mass 

expansion in case of the electromagnetic current and by the 1/m^ term in case of the tensor 
current. The spectra of the correlators of such induced currents start at zero energy. This fact 
makes impossible the standard analysis of the moment sum rules at 0{af) order of perturbation 
theory for n > 4 in case of the electromagnetic current and for n > 3 in case of the tensor 
current. 
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